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1.  STATEMENT  OF  THE  PROBLEM 


In  a  set  of  n  oil  wells  located  in  reservoirs  without  communications,  we 
wish  to  schedule  the  flowrate  of  each  individual  well  in  such  a  way  that  we 
satisfy  the  following  conditions:  cl 

--  the  pressure  drop  of  each  individual  well  is  less  than  or  equal  to  a 
given  value, 

--  the  production  of  any  well  is  bounded  by  lower  and  upper  technical 
limits , 

--  the  total  production  of  the  fields  has  to  satisfy  a  certain  demand 
with  some  specific  technical  limiting  capacities, 

--  the  total  discounted  profit  of  the  oil  production  is  maximized. 


Following  the  basic  fluid  and  cash  flow  mechanisms  of  the  system,  two 
flowrate  arrangements  will  be  considered,  namely: 

—  discrete  flowrate, 

--  continuous  flowrate. 

They  will  be  applied  to  systems  including  single  or  multiwells  through  a 
linear  programming  model  and  an  attempt  to  use  the  Neyman-Pearson  lemma. 


1.1  FLUID  FLOW 


Besides  limitations  on  the  production  of  each  part  of  a  given  oil  field, 
a  special  consideration  has  to  be  made  for  the  maximum  pressure  drop  in  the 
reservoir.  Since  the  maximum  is  located  at  the  bottom  of  the  well  bore, 
this  region  will  have  to  satisfy  at  any  time  the  pressure  constraints. 

The  concept  of  a  well  function  defined  as  the  bottom  hole  pressure  drop 
at  time  t  per  unit  of  constant  flowrate  is  written  as 


Fj(tx)  where  F.(0)  =  0 


1.101 


and  due  to  the  linearity  of  the  diffusivity  equation,  we  can  use  the  super¬ 
position  principle  or  convolution  method  as: 


n 

p0,j  -px,j  -  l  C«i,j 


1.102 


We  want  to  satisfy  for  the  well  j, 


Pn  ,  -  P  ,  Sd, 
0,J  x,j  j 


1.103 


therefore,  in  the  case  of  a  discrete  flowrate, 


ji 

Z( q.  .  -  q,  .  .  )F ,(t  -  t.  ,)  S  6 

i,J  Hi-l,j'  y  x  i,j' 


1.104 


and  without  loss  of  generality,  we  can  assume 


t.  .  -  t.  .  =  At  =  — 

i-l»j  n 


V  y  AtF.(t  -  t.  .)  ?*  d. 

At  .1  x  i,j'  j 


1.105 


1 
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where 

9i  /(t)  =  Um  ftiJ j 

7  n  ->  oo  / 

or  At  -»  0 

The  three  well  functions  which  will  be  used  in  the  calculation  are  repre¬ 
sented  in  figs.  1  and  2. 


Maximum  and  minimum  limitations 

Uniquely  referring  to  the  reservoir  capacities  with  some  pressure  drop 
imposed,  we  can  define  two  limitations  on  the  flowrate,  namely: 

max  constant  flowrate  during  any  time  period  which  will  produce 

the  maximum  permissible  drop  at  the  end  of  the  period. 

max  Qtnin:  constant  flowrate  during  the  entire  time  which  will  produce 

the  maximum  permissible  pressure  drop  at  the  end  of  the  entire 
period. 

For  a  well  j : 

d.  <5 

max  Vx  =  T^t:  “d  Vn  *  FTtT  1.107 

Fj(n  J 


For  a  multiwell 

max  Q  = 
Tnax 


system  which  includes  m  wells: 


1.108 


The  technical  capacities  and  the  demands  must  be  such  that: 

min  (max  q^  ,  min  technical  capacities  well  j)  S  individual  flowrate  ^  1.109 
min  (max  q^^,  max  technical  capacities  well  j)  £  individual  flowrate  qx^  1.110 
min  (max  Q^n,  min  total  technical  capacities)  ^  total  flowrate  1.111 


min  (max  Q^fly>  max  total  technical  capacities  £  total  flowrate  1.112 

and 

demand  S  total  flowrate  1.113 


1.2  CASH  FLOW 


Between  times  t.  and  t.^.  when  the  flowrate  is  q.  .  the  profit  is  such 
i  i+1  K 


that: 


Profit  = 


**  =  f*  1+1  e"XTT  j(t)dt 

tj. 


1.201 


where 

_XTt 

the  discount  function  is  e 

and  where  r  = 

Let 

a  =  AT  and 

-a 

Y  =  e 

or 

t  1  0^  e-XTT 

1 

and 

T  = 

*  Y 

t 

x 

T 


Four  cases  will  be  considered  here: 


Y1 

Y2 

y3 


1.00 

a 

0.90 

0 

0.75 

a 

0.50 

a 

0 

1 

2 

3 


0.000 

0.105 

0.288 

0.693 


The  corresponding  discount  functions  are  shown  on  fig.  3* 


Total  discounted  profit 

a)  Discrete  flowrate : 

t.  and  t .  . : 
l  l+l 


assuming  a  constant  flowrate  q.  .(t)  between 


i+1  "  T 
e 

i 


Profit 


q  (t)dt 
1  )  J 


1.202 
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n 


Total  discounted  profit  =  ^  ^(t)^e  ^  -  e 


i=l 

naT 


1  .  +  T  .  , 

Let  T  =  nAt  and  t  .*  =  — - — 

i  2 


then 


b)  Continuous  flowrate: 


1.203 


1.204 


1.205 

The  profit  function  has  the  dimension  of  dollars  and  it  represents  the  total 
profit  over  the  entire  lifetime  of  the  project. 


2.  INDIVIDUAL  WELL  PROBLEM 

We  consider  one  well  in  a  reservoir  and  it  is  required  that  the  oil 

production  of  the  given  well  be  scheduled  during  a  time  T.  In  a  first  linear 

program  model,  the  total  time  T  is  divided  into  a  series  of  n  time  periods  At  and 

q  .(t)  is  computed.  Then,  the  Neyman -Pearson  lemma  will  be  applied  in  the 
i>  J 

case  of  a  continuous  flowrate  for  the  determination  of  the  exact  function  q j ( t ) . 
It  is  obvious  that,  in  both  cases,  the  cumulative  production,  namely, 


n 


cannot  be  greater  than  the  reserves  of  the  reservoir.  However,  this 


last  constraint  will  not  be  considered  here,  and  we  assume  that  for  any  well 
j  we  satisfy: 


n 


i=l 


or 


f  q.(t)dt  £  R 
JO  J  J 


2.001 


where  R.  are  the  reserves  for  the  well 
J 


j- 


2.1  DISCRETE  FLOWRATE 


From  the  preceding  statements,  in  the  case  of  a  discrete  flowrate  the 

function  q.  .(t)  is  the  solution  of  the  following  system: 

1  >J 


n 


l)  maximize  —  V  e  01 i  q  .(t) 
n  l—i  1  jJ 


i=l 


2.101 


over  all  functions  q.  ,(t)  subject  to 

i  j  j 

2)  the  pressure  constraints: 
n 


or 


Z*i,j(t)W(N  ■  i+iH  -  fj[(n  - i)4t]} 

i=l 


£  <5 


2.102 


1,  ...»  n 


■  i+iH  ■  “j[(n  ■ i)At] } s 

1=1 


2.103 


,  A  2«hk  A 

where  a  .  =  — - —  0  . 

j  n  j 


3)  bounded  by  some  technical  limiting  capacities  as: 


q.  .min(t)  £  q.  .(t)  £  q.  .max(t) 
i ,  j 


2.104 
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4)  where  the  entire  lifetime  of  the  project  is  defined  as: 

0  5  tL  §  T  2.105 

This  problem  is  easily  solved  by  the  use  of  computers  and  well-improved 
lineal'  programming  techniques.  The  solutions  obtained,  however,  show  a  de¬ 
pendence  under  the  selected  value  of  n.  The  true  maximum  could  be  obtained 
by  greatly  increasing  n.  But  the  size  of  our  problem  is  fast  becoming  large 
and  lengthy.  It  is  therefore  interesting  to  try  to  determine  the  actual  true 
maximum  by  another  method  dealing  with  continuous  flowrates. 


0  5  t  £  T 


2.204 


10 


The  above  system  will  be  somewhat  modified  through  a  change  of  vari¬ 
ables,  in  order  to  take  advantage  of  the  Neyman-Pearson  lemma  as  presented 
in  appendix  5.1  and  by  fig.  14. 


Let 

t 

T  T 

then 

0  *  x  *  1 

2.205 

and 

t) 

q.(x)  =  - — 

J  q.max 

J 

2.206 

or 

0  £  q.(x )  S  1 

2.207 

Now  consider  the  equation  2.201  which  is  equivalent  to 

rT  “I 

minimize  /  -  e  q.(t)dt 

J  0  J 


Or,  if  we  let 


a(x)  =  -e 


-ax 


2.208 


the  objective  function  of  the  system  is 


minimize  /  a(x)q.(x)dx 
J  n  J 


2.209 


Consider  the  function  l(t)  defined  as 


I(t)  =  [  q.'(t)u>  (T  -  t)dt 
J  0  V  V 


and  let  us  transform  this  function  in  the  following  way: 


Mt)  =/0V(TT)“d[T( 


1  -  x) 


Tdx 


dq.(x) 

—A _ 

dx 


.  dq  (x)  q.max 

q,  '(TTjt^lj'Crr)  =  -k - 


q  .max 

j 
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and,  finally: 


f  T  dq  (t)  _ 

I,  (x)  =  q.max/  — —  w 
1  j  dx 


T(  1  -  t) 


dT 


2.210 


Integrating  by  parts,  equation  2.202  can  be  stated  as 


VT)“j 


T(  1  -  t) 


T  pT  dW 

-  q(T) 

0  ^0 


4 


T(  1  -  T) 


ST 


A  . 

dr  5  — L 


q.max 

J 


2.211 


Let 


dw. 

b(x)  =  -  -i 


T(1  -  t) 


dx 


2.212 


and 


c(t)  '^■{qj(T)“ifT,1'T)]l'  !o 


2.213 


or 


c(t)  =  ii  +  'ij<°)“i(T)-qj(T)“j[Tfl  -T)  j 


2.214 


Our  initial  system  is  now  the  following: 


minimize  /  a(x)q.(x)dx 

J0  3 


2.215 


over  all  functions  q.(t)  subject  to  the  constraints: 

J 


a)  0  S  q  (x)  S  1,  0  5  t  5  1 
J 


2.216 


b)  [  b(x)q.(x)dx  £  c(x) 


0 


2.217 
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where  a(x),  b(T),  and  c(x  )  are  given  functions  such  as 


a(x)  =  -e 


-C(X 


dw . 

b(x)  =  -  -J 


T(1  -  t) 


dx 


A. 

c(x  )  =  — - —  +  q.(0)ui.(T)  -  q.(x)w. 
q  max  V  J  V  '  J 

J 


T(  1  -  t) 


2.218 

2.219 

2.220 


This  form  of  the  system  allows  free  use  of  the  Ileyraan -Pear son  lemma, 
from  which  we  obtain  the  optimum  solution. 


2.3  PROPOSED  SOLUTIONS 

A  numerical  solution  is  now  presented  for  the  case  of  an  individual  well 
surrounded  by  a  radial  circular  reservoir.  The  derivation  of  the  well  func¬ 
tion  of  such  a  system  is  shown  in  appendix  5*2.  Three  different  wells  will 
be  considered.  In  this  first  part  only  one  single  well,  say  well  j,  is  in¬ 
cluded  in  the  system  that  we  want  to  schedule. 

All  the  units  in  this  study  are  expressed  in  terms  of  the  centimeter- 
gram-second  system,  the  results  being  shown  in  practical  units. 


Characteristics  of  the  wells 

3 

1)  thickness  of  the  formation:  h  =  10  cm 

2)  well  radius:  a  *  10  cm 

3)  outside  reservoir  radius:  R  =  2  y  lO'*  cm  or  R  5  1  mile 

4)  porosity:  u>  =  10$ 

5)  permeability: 

well  (1):  k  =  10-12  CGS  or  0.1  millidarcy 

well  (2):  k  =  10-11  CGS  or  1.0  millidarcy 

well  (3):  k  -  10~10  CGS  or  .10.0  millidarcy 
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10  2 

6)  compressibility  of  the  fluids:  p  =  10  CGS  (cm  per  dyne) 

7)  viscosity  of  the  fluids:  ji  =  10  2  CGS  (Poise  or  dyne-second 

2  v 

per  cm  ) 

g 

The  maximum  pressure  drop  allowed  here  is  of  10  dynes  or  1,450  psi. 

The  numerical  solution  presented  covers  a  time  interval  of  1  year.  As 
can  be  seen,  the  flowrate  profiles  are  independent  of  the  economic  conditions 

imposed  in  the  system. 

2-.  31  DISCRETE  FLOWRATE  LINEAR  PROGRAMMING  MODEL 

We  want  to  solve  equations  2.101,  2.102,  and  2.103  for  a  well  j,  where 

T  =  1  year  and  n  =  12 


Well  function 

In  this  case  the  function  u^{e)  (see  appendix  5.2,  equation  5.211)  can 
be  used: 

3  (e)  =  |(Log  e  +  0.80907) 


Or  let 


F.(t)  =  a  log  (t)  +  b 

u 


2.3101 


where : 


a 

'B 

well  (l) 

1.83239  x  106 

-  1.18852  X  10fc 

well  (2) 

1.83239  x  105 

0  6438  x  I0h 

well  (3) 

1.83239  X  104 

2.4762  x  10* 
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Maximum  and  minimum  flowrates 


From  equation  1.107  and  when 


6. 

J 


10  dynes : 


F.(T) 

J 

Fj(t/6) 

F.(T  '12) 

J 

Maximum 
(CGS ) 

Max  a 

Tnax 

EEU 

BBS 

n  =  12 

well  (1) 

0EHE3 

7.974 

7,974 

8.997 

9.466 

69.575 

69.575 

77.238 

’  0 

well  (3) 

637.055 

617,055 

676 . 589 

Assuming  that  the  technical  capacities  are  such  that 


minimum  technical  capacities  =  0.0 


and 


maximum  technical  capacities  >  max  q.max 


equations  1.109  and  1.110  are  reduced  to 


0  < 


flowrate  q.  5  max  q.max 
J  .] 


2.3102 


Solution 

The  solution  of  this  simple  linear  programming  problem  is  easily  obtained 
by  using  the  Share  program  number  H0SCM3  on  the  7094  IBM  computer. 

In  fig.  4  a  case  is  presented  where  upper  and  lower  limits  exist  on  the 
individual  flowrates  for  well  (2),  assuming  four  different  discount  functions, 
namely,  aQ,  o^,  a^,  a^,  as  defined  in  section  1.2. 

The  optimum  flowrate  expressed  as  the  ratio  q  average  q  max  for  differ¬ 
ent  values  of  q  max  'max/a  is  shown  on  fig.  5  for  the  three  wells.  The 

'  d'lax 

profit  is  also  represented,  and  the  influence  of  the  number  n  on  the  final 


<  * 


results  can  be  seen. 
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2.32  CONTINUOUS  FLOWRATE  NEYMAN -PEARSON  IEMMA 


From  equations  5.201  and  2.219, 


b(x)  =2 


n=l 


P  2a2 

n 


T(  1  -  t  ) 


2.3201 


from  equation  2.220, 


c(T)  =  53^  +  «3<°)[IO« R/  -  «  Zv 


P  2a2 
n 

K 


T*l 


r 


-  qj(T)  Log  R'  -  2  ^«ne 


P  2a2 
rn 

K 


T(1  -  T  )*l 


2.3202 


and  from  equation  5.1104, 


A(t) 


■  i  X  V 


„  2  2 
0  a 
rn 

K  T(1  -  t  ) 


-Iv 


ft  2  2 

Pn  a 
K 


2.3203 


These  three  functions,  b(x),  c(t),  and  A(t  )  are  well  defined  for  a  given 
well  and  present  some  interesting  properties.  For  instance,  b(r)  being  a 
monotonically  non-negative  increasing  function  with  a  finite  value  when 
t  =  0,  and  a(x)  having  also  a  finite  value  when  t  =  0,  the  supremum  value 
of  cXq',  as  defined  in  appendix  5«1>  is  such  that 

°o/  =  ■  hi 


2.3204 
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Therefore,  x^  the  intersection  of  a(x)  and  a  'b(x),  is  equal  to  one. 

Thus  a  given  well  will  be  flowing  at  its  maximum  flowrate  from  t  =  0 
to  x  =  T2*  T2  ^efined  as  the  intersection  of  A(x)  and  c(x).  Then  it 

is  only  necessary  to  determine  x^  and  to  verify  that  A(x )  £  c(x)  when 

0  s  T  s  v 

A. 

If  v:e  notice  also  that  A(0)  =  0  and  c(0)  =  - —  >  q  then 

q  .max 
J 

q^O)  =  1  2.3205 

This  means  that  the  well  has  to  be  set  at  its  maximum  value  at  zero  time. 

This  maximum  value  will  remain  up  to  the  time  x^,  and  the  intersection  of 

A(x)  and  c(x)  can  be  determined  by  setting  q.(x)  =  1  in  equat'on  2.3202. 

J 

However,  a  careful  observation  of  the  problem  shows  that  we  have  to  satisfy 
equation  2.217,  where  one  of  the  limits  is  a  function  of  x.  Therefore,  in 
some  cases,  it  could  be  hazardous  to  apply  the  Neyman-Pearson  lemma  system¬ 
atically.  A  more  complete  method  of  optimal  control,  as  the  maximum  princi¬ 
ple  of  Pontriagine,  would  be  preferable  (ref.  5,6).  This  limitation  restricts, 
to  a  certain  extent,  a  simple  application  of  the  Neyman-Pearson  lemma  for  this 
particular  problem. 

3  MULTIWELL  SYSTEM 

The  multiwell  system  is  considered  here  as  an  extension  of  the  individual 
case  when  m  wells  are  producing  from  a  given  set  of  reservoirs.  When  more 
than  one  well  acts  on  the  same  system,  however,  interference  effects  are  fre¬ 
quently  observed.  These  effects  will  not  be  taken  into  account  here;  only 
the  simplified  case  where  all  the  wells  are  independent  will  be  developed, 
assuming  that  each  well  has  the  same  effect  as  if  it  were  alone  in  the  system. 
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Therefore,  the  cumulative  production  will  consist  of  the  linear  summation  of 
all  the  wells  and  the  total  cash  flow  will  mediate  the  corresponding  produc¬ 
tion  with  its  economics. 


3.1  DISCRETE  FLOWRATE 


Initially,  it  was  thought  that  a  generalization  of  the  fundamental  lemma 
of  Neyman  and  Pearson  might  be  used  for  the  treatment  of  the  multiwell  system 
(ref.  4).  However,  due  to  the  above-mentioned  limitations  regarding  the  indi¬ 
vidual  well,  this  method  had  to  be  used  with  caution.  It  would  be  more  advis¬ 
able  to  use  the  optimum  control  principles  presented  in  references  5  or  6. 

Thus,  this  study  will  be  concerned  only  with  some  aspects  of  the  discrete 
case  where  the  problem  to  be  solved  is  the  following: 

--  Find  the  optimum  flowrate  of  a  multiwell  system  including  m  wells, 
with  the  following  constraints: 

--  pressure  drop  maximum  at  every  well 
—  upper  and  lower  bounds  on  the  wells 
--  a  certain  demand  to  be  satisfied 
Then  the  formulation  of  the  problem  can  be  set  as: 
l)  maximize 

'  m  n  „ 

l  n  ^Ti 

0=1  i=l 

over  all  functions  q.  .(t)  subject  to 

1  >  J 


3.1 


m  n 

2)  L  L  qi,J(t)  {SJ 

j=l  i=l 


(N  -  i  +  l)At 


-  u> . 

JL 


(N  -  i)At 


H, 


N  =  1,  . . . ,  n 

3)  bounded  by  some  technical  limiting  capacities  as 


3.2 


q.  .min(t)  S  q,  .(t)  S  q  max(t) 
1  >  J  1  *  J  1  >  J 


3.3 
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f 

m 

U)  I  li.jC-)  s  «*) 

J-l 

5)  where  the  entire  lifetime  of  the  project  is  defined  as 

0  £  t  S  T 

3.2  PROPOSED  SOLUTION:  LINEAR  PROGRAMMING  MODEL 

This  linear  programming  model  has  been  solved  for  different  cases. 

The  basic  matrix  is  presented  in  fig.  6,  where: 

n  =  6,  m  =  3j  q,  .min  =  0 

y  J 

Q 

q.  .max  =  »,  Q( t )  =  constant,  A  .  =  10  dynes 
1  >  J  J 

Due  to  the  number  of  constraints,  the  size  of  the  matrix  increases  as 
n  is  increased.  So,  when 

n  =  1  we  get  a  matrix  of  8  rows  x  l6  columns 
n  =  6  we  get  a  matrix  of  28  rows  x  46  columns 
n  =  12  we  get  a  matrix  of  52  rows  x  82  columns 
and  when  n  is  large,  the  problem  is  greatly  simplified  by  using  the  decomposi 
tion  principle  for  linear  programming. 

Also,  the  number  of  rows  and  columns  in  the  basic  matrix  determines 
whether  the  primal  or  the  dual  type  of  problem  should  be  used,  the  fewer 
the  number  of  rows  in  the  matrix,  the  easier  to  solve  the  problem. 

An  obvious  simplification  can  be  obtained  by  fixing  a  priori  n  *  1. 

As  can  be  seen,  however,  a  true  maximum  is  not  attained  by  this  method. 

But  due  to  the  small  increase  of  the  total  cash  flow  (a  few  per  cent)  when 
using  n  =  6  instead  of  n  =  1,  this  simplification  can  be  worthwhile.  It 


3.4 
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can  also  be  noted  that  by  changing  the  value  of  n,  it  is  possible  to  modify 
the  well  recovery  to  a  considerable  extent.  This  is  illustrated  by  fig.  8, 
case  B^»  and  fig.  13,  case  B 

Two  categories  of  results  are  shown  from  figs.  7  to  13. 
l)  Constant  output-equality  cases  where 


m 

l 

j=l 


q.  .(t)  =  constant  (in  arbitrary  units) 

1  9  J 


3.6 


a)  no  limitations  on  individual  flowrate 
fig.  7  n  ■  6 

fig.  12  n  -  12 

b)  limitations  on  individual  flowrate 
fig.  11  n  =  6 

2)  Total  bounded  output -inequality  cases  where 


m 


I 


q.  .(t)  5  constant  (in  arbitrary  units) 
1 » J 


3.7 


J=1 

figs.  8,  9,  10  n  =  6 
fig.  13  n  =  12 

The  results  are  shown  with  the  following  characteristics: 

--  on  the  left  side  of  the  page,  the  economics  based  on  the  a-coeffi- 
cient  values 

—  in  the  center  of  the  page,  the  production  profile 

—  on  the  right  side  of  the  page,  the  total  and  individual  output  for 
the  three  wells 

Four  basic  mechanisms  have  been  observed: 

A.  For  the  constant  output-equality  cases,  an  improvement  of  the  econom¬ 
ics  leads  to  an  increase  of  the  cash  flow  associated  with  a  decrease  of  the 


recovery  from  each  well. 
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B.  For  the  total  bounded,  output -inequality  cases  we  fix  a  constant 
demand  and  follow  the  change  of  the  optimum  flowrates  for  different  values 
of  the  constant.  The  interesting  parameters  are  the  decrease  of  the  demand 
and  the  change  of  the  economics. 

Then: 

--  for  large  values  of  the  demand,  all  the  wells  have  to  flow  to 
their  maximum; 

--  when  decreasing  the  demand,  the  best  wells  are  kept  at  their 
maximum  value  while  the  other  wells  are  reduced  conversely  proportional  to 
their  productivity; 

--  for  a  given  value  of  the  demand  corresponding  to  a  maximum 
constant  demand  for  the  whole  period,  the  situation  is  completely  modified; 

--as  soon  as  the  demand  is  below  the  previous  value,  the  bad 
wells  are  fixed  to  their  maxima  and  the  final  fit  is  made  on  the  best  well 
in  reference  to  its  productivity. 

This  procedure  can  be  followed  on  figs.  8,  9,  10,  and  13-  On  figs.  7  and 
12,  the  decrease  of  the  maximum  value  of  a  constant  demand  is  shown  to  be  a 
function  of  the  different  economics  imposed  on  the  wells. 

C.  Effects  of  upper  and  lower  bounds  are  shown  on  fig.  11,  which  can  be 
compared  with  fig.  7- 

D.  Production  profile.  During  the  entire  period  the  well  production  has 
to  satisfy  a  certain  profile .  The  first  steps  of  this  profile  are  especially 
interesting. 

The  results  shown  can  be  considered  as  typical  behavior  for  a  three-well 
system  where 

well  (l)  is  a  low  productivity  well 
well  (2)  is  a  medium  productivity  well 
well  (3)  is  a  high  productivity  well 
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The  different  values  of  the  economics  must  be  sufficient  to  cover  most  of 
the  actual  cases . 

Using  an  IBM  709^  computer  and  Share  program  number  H0SCM3,  the  com¬ 
puting  time  was  2  seconds  for  each  case  when  n  =  6. 

4  CONCLUSION 

--  A  good  estimate  of  the  optimum  cash  flow  can  be  obtained  with  r  =  10. 

--  The  optimum  recovery  of  a  well  is  very  sensitive  to  the  economics  im¬ 
posed  on  the  system,  as,  for  example,  the  lower  and  upper  bounds  restrictions. 

—  The  bad  wells  must  produce  at  a  flowrate  proportionately  higher  than 
the  best  wells. 

--  The  production  profile  of  the  well  is  not  obvious  in  most  cases,  but 
an  accurate  approach  to  the  optimum  profile  can  be  obtained  at  moderate  cost 
by  using  a  linear  programming  model. 
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5  APPENDIX 

5.1  THE  NEYMAN -PEARSON  LEMMA 


5.11  Statement  of  the  problem 

We  wish  to  minimize  a  linear  functional  of  the  form: 


over  all  functions  f(x  )  subject  to  the  following  constraints: 

a)  0  §  f(x)  *  1,  0<t<1 

b)  [  f(x)b(x)dx  5  c 


5.1101 


5.1102 

5.1103 


where  a(x)  ^  0  and  b(x)  ^  0  are  given  functions,  and  c  ^  0  is  a  known  constant. 


5.12  Solution 

The  solution  is  determined  as  follows.  Let  set  functions 
E'  =  E~(a'),  E  =  E(a'),  E+  =  E+(a') 


be  defined  for  <  a  <  +»  as : 

E"*(a')  =  [x;  a(x)  <  a'b(x)] 

E(a')  *  [x;  a(x)  =  a'b(x)] 

E+(a/)  =  [x;  a(x)  >  a'b(x)] 

Determine  by  the  condition  that  be  the  supremum  over  all  non¬ 
positive  a'  satisfying  the  inequality 

f  b(x)dx  5  c 
JE~ 


5 .  HOU 
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and  let 


E'(o0')  =  E0‘ 
E  (“o')  ■  E0 

EV>  -  V 


Then  the  set  of  minimizing  functions  f  is  given  by: 

a)  f  (t)  =  1  on  Eq" 

b)  f  (t)  =  0  on  Eq+ 

c)  f  (t)  =  arbitrary  on  EQ,  satisfying  only  the  condition 


0  <  f  (x)  <  1,  0  *  x  5  1 


and 


T  1  * 

I  f  (t  )b(x  )dx 


0 


c  if  o  '  <  o 


or  the  conditions 


0«f  (t)?1,  0  ?  x  *  1  if  a  '  =  0 


5 . 13  Outline  of  the  procedure 

a)  choose  a '  =  a . 1  <  0 

i 

b)  find  E~(ai/)  =  [x;  a(x)  <  o^b(x)] 

c)  compute  A  ■L  b(x  )dx 

E" 

d)  if  A^  >  c,  decrease  a^'  and  go  to  (a) 

if  A/  \  <  c,  increase  a.'  and  go  to  (a) 

[T  )  i 

Ontimal  is  such  that  for  +  Ao^,  A^  >  c, 


See  fig.  lU. 
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5.2  WELL  FUNCTION  OF  THF  SINGLE  WELL 


The  well  function,  defined  as  the  pressure  drop  in  a  given  well  corre¬ 
sponding  to  a  constant  unit  flowrate  at  time  t,  can  be  computed  in  many  ways 
(see,  for  instance,  ref.  2).  Here  we  are  interested  in  the  analytical  formu¬ 
lation  of  the  well  function  corresponding  to  a  radial  system,  and  we  show  that 
such  a  function  can  be  written  as: 


5.201 


where  each  term  is  defined  below. 


Radial  flow  of  a  tomogeneous  fluid  into  a  single  well 

The  well  completely  penetrates  the  reservoir  and  the  fluid  flows  uni¬ 
formly  in  all  directions  radiating  to  the  axis  of  the  well  bore.  The  system 
has  axial  symmetry  and  we  can  describe  the  fluid  pressure  variations  by  the 
Laplace  equation  in  cylindrical  coordinates,  namely: 

1  _  _1  c3p  5.202 

~2  +  r  5r  "  K  5r 
or 


with  dimensionless  variables  as  follows: 


distances:  r'  *  — 

a 

time:  6  =  — 

a 

pressure:  w  =  (P 


P) 


2nhk 

HQ 


Equation  5.202  is  then  reduced  to: 


5^  to  1  du) 
dr'2  "  ?75F7 


du) 

56 


5.203 


Initial  conditions 

—  constant  pressure  in  the  whole  reservoir  at  the  start,  or: 

6=0,  w  =  0  5.204 


--  boundary  conditions: 

constant  flowrate  from  the  well,  or: 

)  ■  -1  for  any  6  5.205 

r'=l 

no  flowrate  from  the  outside  border  and  constant  pressure  at  the 
reservoir  limit,  or: 

u  =  0  for  any  6  5.206 

Analytical  solution 

The  solution  of  the  above  system  with  the  given  initial  and  boundary 
conditions  is  such  that: 


where  the  B  's  are  roots  of: 
rn 
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If  we  notice  that  the  last  terra  of  5.207  decreases  rapidly  as  the  time  in¬ 
creases,  after  a  certain  time: 

<^(6)  =  Log  R 7  5.209 

Equation  5-207  represents  the  unsteady  state  and  equation  5.208  represents 
the  steady  state.  <3  (6)  is  shown  on  fig.  15. 

Approximations  w ^(6),  w  (6 )  and  ) 

On  fig.  15  we  can  differentiate  three  zones,  namely: 
zone  1:  the  initial  part  of  the  curve, 

zone  2:  the  central  part  presenting  a  straight  semilog  variation, 
zone  3:  the  end  portion  with  the  asymptotic  value  Log  R ' . 

These  three  zones  are  bounded  by  times  6^  and  6  which  will  be  considered 
further. 

If  we  now  vary  R',  we  see  on  fig.  16  that  all  the  curves  representing 
are  tangent  to  the  curve  for  which  the  external  radius  is  infinite. 

Physically  this  property  is  very  easy  to  justify.  Then  the  calculus  shows 
that  when  the  radius  is  infinite: 

5.210 

and  for  values  of  6  larger  than  6 this  last  expression  converges  to 

5.211 

The  semilogarithmic  variation  for  the  case  of  a  reservoir  of  finite 
radius  can  be  used  up  to  a  certain  value  B^  which  will  be  considered  further. 
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If  we  now  assume  that  the  well  radius  is  null  and  the  outside  radius 
infinite,  the  solution  is  such  that: 


<*>4(0) 


5.212 


This  last  function  can  be  of  some  use  when  we  want  to  conroute  the  well  func¬ 
tion  in  an  interval  such  that  6  5  6 .  Besides  its  obvious  simplicity  compared 
with  )  or  w  (0),  an  interesting  use  of  the  Wq(6)  function  can  be  seen  through 
its  derivative  when  referring  to  the  Neyman -Pearson  method,  for  instance. 


Interval  for  the  use  of  6^  and  6 ^ 

6  __  — 

1:  the  convergence  of  u>^(0)  to  w  (0)  is  illustrated  on  fig.  17,  where 

we  can  see  that  the  approximation  w^(6)  for  u>^(6)  is  excellent  as 

2  / 

far  as  6  §  10  and  this  limit  is  practically  independent  of  R  . 

6 

2:  the  value  of  6^  is  defined  as  the  limiting  value  of  6  ut>  to  which 
we  get  an  appropriate  convergence  between  u^(0)  and  u>^(0).  If  we 
fix  a  ratio  of  convergence  such  that: 

w,(e)  -  u  (e) 

1  _  ,rr -  »  6/1000  5.213 

we  get  fig.  18,  from  where: 

Log  6^  =  8.  Log  R/  +  b  5.214 

which  shows  that  C  is  dependent  on  R/  under  a  simple  formulation. 
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Calculation  of  a  and  i 
n  '  n 


c*n:  from  equation  5.207, 


J0  <^nR/> 

f,n2|Jl2(^n>  -J02"!nR') 


>,2(:U 


>  2|  1  v  n 

n  L.T  2f  n  R / 


-To<V'> 


5.215 


but  from  5. 208, 


YAv 

V02(flnR') 


5.216 


and  for  small  values  of  p 


Log  =  -  Log  p  * 


(A  ■  A  s 


when  R'  >  100  ^  1 


Y  is  large,  say  >  600 


Y  (PR')  is  small,  say  <  1 


<*n  =  2.46738  Y^^R') 


5.218 


29 


When  R'  >  100,  p  is  small  and  J.(0)  =  0;  then  J..(p  )  —  0,  and 
n  1'  ’  In 

5.222 

As  soon  as  n  is  larger  than  15,  the  positive  roots  of  5.222  are  given  annroxi- 
mately  by: 

n  £  15  5.223 

The  following  table  summarizes  the  values  of  and  p^R'  -  100 )  obtained 
from  the  above  considerations. 
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Table  1 


n 

a 

n 

(3  R' 
rn 

n 

0 

n 

P  R' 

'  n 

n 

a 

n 

P  R' 
rn 

1 

0.6415 

2.405 

11 

0.0465 

33.76 

21 

0.0241 

65.18 

2 

0.2834 

5.520 

12 

0.0425 

36.92 

22 

0.0230 

68.33 

3 

0.1812 

8.654 

13 

0.0392 

40.06 

23 

0.0220 

71.47 

4 

0.1331 

11.79 

l4 

0.0364 

43.20 

24 

0.0211 

74.62 

5 

0.1051 

14.93 

15 

0.0339 

46.34 

25 

0.0202 

77.76 

6 

0.0869 

18.07 

16 

0.0317 

49.48 

26 

0.0194 

80.90 

7 

0.0740 

21.21 

17 

0.0299 

52.62 

27 

0.0187 

84.04 

8 

0.0645 

24.35 

18 

0.0281 

55.77 

28 

0.C180 

87.18 

9 

0.0571 

27.50 

19 

0.0267 

58.90 

29 

0.0174 

90.32 

10 

0.0512 

30.63 

20 

0.0253 

62.05 

30 

0.0168 

93.46 

The  asymptotic  well  function  u>q(  B ) 


In  some  particular  cases,  and  for  a  specific  interval  of  time,  the  use 
Of  u>2(e),  w^(0),  or  can  simplify  the  problem.  This  is  especially  ob¬ 
served  when  using  the  function  between  6^  and  6^.  This  last  interval 

is  sufficiently  large  in  actual  so  that  a  semilogarithm  variation  is 

justified  almost  a  priori  for  a  given  well.  Nevertheless,  it  seems  interest¬ 
ing  to  refer  to  a  well  function  covering  the  complete  time  scale.  Such  a 
function  can  be  estimated  from  the  computed  coefficients  and  P^.  Only 
finite  values  of  n  will  be  used  in  order  to  simplify  the  formulation  without 


altering  its  effective  accuracy.  Due  to  its  property  this  function  can  be 
named  the  asymptotic  well  function  w^(G). 
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A  computation  of  u5Q(e)  using  the  a n  and  Pn  coefficients  from  table  1 
shows  that  the  results  reproduce  those  presented  in  reference  7  and  that 
for  R'  >  100,  the  convergence  of  the  number  of  terms  is  the  following: 

6  =  1  n  =  8l 

6  =  10  n  =  29 

e  =  102  n  =  10 

6  =  103  n  =  U 

with  a  rapid  convergence  to  one  for  larger  values  of  6  depending  on  tne 
actual  value  of  R/. 

It  is  then  suggested  that  only  the  first  four  values  of  n  be  introduced 
for  obtaining  a  good  estimation  of  the  well  function  of  a  single  well  located 
at  the  center  of  a  circular  reservoir. 


NOMENCLATURE 


a(x) 

b 

b(x ) 
c(x) 
Et(-t) 
fj(T) 

'!<*)■ 

h 

i 

I-l(x) 

J 

k 

K: 


m 

n 


N 


0,3 

Pi,J 

pj(t) 


x,j 


*±,i 

qj(t) 

q/(t) 

q^max 

q  .min 

J 


well  radius 
function  (see  2.208) 
coefficient 
function  (see  2.212) 
function  (see  2.213) 
exponential  integral 

modified  pressure  drop  at  well  j  for  modified  time  x 

pressure  drop  at  well  j  for  a  constant  flowrate  during  time  t 

reservoir  thickness 

integer  =  1,  2,  . ..,  n,  n  +  1 

function  (see  2.210) 

integer  =  1,  2,  ...,  m  -  1,  m 

reservoir  permeability 

hydraulic  diffusivity 

number  of  wells  in  the  multiwell  system 

number  of  equal  intervals  in  period  T 

integer  =  1,  ...,  n 

initial  pressure  at  well  j 

pressure  at  well  j,  time  t^ 

pressure  at  well  j,  time  t 

pressure  at  well  j,  time  x 

discrete  flowrate  at  well  j,  time  t^ 

continuous  flowrate  at  well  j,  time  t 

derivative  of  flowrate  at  well  j,  time  t 

maximum  flowrate  allowed  on  well  j 

minimum  flowrate  allowed  on  well  j 

flowrate  at  well  j,  time  t 
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^>0 

Q(t) 


S' 


m 


R 


R  . 
J 


Y01YU>J0>Ji 


a 

a 

n 

P 

Pn 

Y 

*i 

A. 

J 

\ 

6 

\ 

M 

T 

* 
T  . 

1 

U) 

w(e) 


derivative  of  flowrate  at  well  j.  time  t 

x 

total  flowrate  of  the  reservoir  at  time  t 
total  production  of  the  well  m 

reduced  well  radius 
reservoir  radius 
reserves  for  well  j 
time  0  ^  t  5  T 

total  period  under  consideration 

time  at  which  the  flowrate  starts  0  £  t^  5  T 

time  OSt  ^  T 
x 

integer  =  1,  2,  . . . ,  n 
Bessel  functions 
function 

roots  (see  5*220) 
oil  compressibility- 
roots  (see  5*223) 
coefficient 

actual  pressure  drop  allowed  at  well  j 

modified  pressure  drop  allowed  at  well  j 

interval  of  time 

reduced  time 

coefficient 

oil  viscosity 

modified  time 

average  modified  time 

reservoir  porosity 

well  function 
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